Pressure dependence of the Curie temperature in Ni 2 MnSn Heusler alloy: A 

first-principles study 
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The pressure dependence of electronic structure, exchange interactions and Curie temperature 
in ferromagnetic Heusler alloy Ni2MnSn has been studied theoretically within the framework of 
the density-functional theory. The calculation of the exchange parameters is based on the frozen- 
magnon approach. The Curie temperature, T c , is calculated within the mean-field approximation 
by solving the matrix equation for a multi-sublattice system. In agrement with experiment the 
Curie temperature increased from 362K at ambient pressure to 396 at 12 GPa. Extending the varia- 
tion of the lattice parameter beyond the range studied experimentally we obtained non-monotonous 
pressure dependence of the Curie temperature and metamagnetic transition. We relate the theo- 
retical dependence of T c on the lattice constant to the corresponding dependence predicted by the 
empirical interaction curve. The Mn-Ni atomic interchange observed experimentally is simulated to 
study its influence on the Curie temperature. 

PACS numbers: 75.50.Cc, 75.30.Et, 71. 15. Mb, 74.62.Fj 
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I. INTRODUCTION 

The pressure dependence of the Curie tempera- 
ture provides important information on a ferromag- 
netic system and is an object of intensive studies 
both experimental 1 ! 2 : 3 ! 4 and theoreticaL 5 ! 6 ! 7 ! 8 ' 9 ! 10 The 
key question here is the character of the variation of var- 
ious magnetic properties with decreasing distances be- 
tween magnetic atoms. In an early work, Castcllitz 11 
proposed an empirical rule (interaction curve) that de- 
scribes the dependence of the Curie temperature of the 
Mn-containing ferromagnetic alloys with 4-5 valence elec- 
trons per molecule on the ratio R/d where R is the 
nearest-neighbor Mn-Mn distance and d is the radius 
of the atomic Mn 3d shell. The curve is supposed to 
represent the Curie temperatures of various systems at 
ambient pressure as well as the pressure dependence of 
T c of a given system. The function is not monotonous 
and has a maximum at the R/d value of about 3.6 (see 
Fig. [SJ. According to the interaction curve, one can 
expect dT c /dP > for alloys with R/d > 3.6 (e.g., 
Ni2MnSn and CuaMnln). On the other hand, the sys- 
tems with R/d < 3.6 (e.g., NiAs-type MnAs, MnSb and 
MnBi) are expected to have negative pressure depen- 
dence of the Curie temperature. These predictions are 
in agreement with experiment 

Recently Kanomata et al. suggested a generalization 
of the interaction curve to the case of 6-7 valence elec- 
trons per chemical formula^ These systems form a new 
branch of the dependence of the Curie temperature on 
the Mn-Mn distance (Fig.©. The available experimen- 
tal values of the pressure derivative of the Curie temper- 
ature, dTc/dP , for Heusler alloys are consistent with 
those expected from the interaction curveiiSiiLiS 

Early experiments on the pressure dependence of the 
Curie temperature of Heusler alloys have been performed 
in a low pressure region (less than 0.5 GPa). Recently 
Gavriliuk et ali^ have studied structural and magnetic 



properties of N^MnSn in the pressure interval up to 10.8 
GPa. The authors have found an increasing linear depen- 
dence of the Curie temperature on applied pressure. The 
Mossbaurer spectroscopy revealed partial interchange of 
the Mn and Ni atoms. 

The purpose of the present work is a first-principles 
study of the electronic structure, exchange interactions 
and Curie temperature in Ni2MnSn as a function of pres- 
sure. The main attention is devoted to the interval of 
the interatomic Mn-Mn distances from 4.26A to 4.06A 
that corresponds to the available experimental variation 
of this parameter. These values of the Mn-Mn distance 
are far above the value of 3.6A that, according to in- 
teraction curve, separates the regions of positive and 
negative pressure gradients of the Curie temperature for 
this group of systems. To verify the appearance of the 
non-monotonous behavior we extended the calculation 
to smaller values of the lattice constant corresponding to 
larger applied pressures. We compare empirical and cal- 
culated dependencies. The influence of the Mn-Ni atomic 
interchange on the magnetism of the system is also stud- 
ied. 

The paper is organized as follows. In Sec. II we present 
the calculational approach. Section III contains the re- 
sults of the calculations and discussion. Section IV gives 
the conclusions. 



II. CALCULATIONAL METHOD 

The calculations are carried out with the augmented 
spherical waves method^ within the atomic-sphere 
approximation^ The exchange-correlation potential is 
chosen in the generalized gradient approximation^ A 
dense Brillouin zone (BZ) sampling 30 x 30 x 30 is used. 
To establish the relation between the lattice parameters 
and applied pressure we use the following expression ob- 
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tained experimentally in Reffl9l 

^M = -aP + 6P* (1) 

where a = 8.64 • lO^GPa' 1 , b = 1.13 ■ lO^GPa" 2 and 
Vo is the volume of the unit cell at the ambient pressure. 
The radii of all atomic spheres are chosen equal. 

We describe the interatomic exchange interactions in 
terms of the classical Heisenberg Hamiltonian 

#e// = -E E J RR S R S R ( 2 ) 

Mi w R,R' 

(pR^r/R') 

In Eq.Q, the indices \x and v number different sublat- 
tices and R and R' are the lattice vectors specifying the 
atoms within sublattices, is the unit vector pointing 
in the direction of the magnetic moment at site (fi, R). 
The systems considered contain three 3d atoms in the 
unit cell with positions shown in Fig^ 

We employ the frozen-magnon approachS&Si^ to cal- 
culate interatomic Heisenberg exchange parameters. The 
calculations involve few steps. In the first step, the ex- 
change parameters between the atoms of a given sublat- 
tice fi are computed. The calculation is based on the 
evaluation of the energy of the frozen-magnon configura- 
tions defined by the following atomic polar and azimuthal 
angles 

°r = > ^ = q-R+0 M - (3) 

The constant phase </) M is always chosen equal to zero. 
The magnetic moments of all other sublattices are kept 
parallel to the z axis. Within the Heisenberg model (J2J) 
the energy of such configuration takes the form^i 

E»» (0, q) = Eg" (6) + sin 2 6J^ (q) (4) 

where Eq m does not depend on q and the Fourier trans- 
form J Miy (q) is defined by 

J^(q)=^J ^exp(zq.R). (5) 

R 

In the case of v = fi the sum in Eq. (JSJ) does not 
include R = 0. Calculating E^^(Q, q) for a regular q- 
mesh in the Brillouin zone of the crystal and performing 
back Fourier transformation one gets exchange parame- 
ters for sublattice [i. 

The determination of the exchange interactions be- 
tween the atoms of two different sublattices \i and v is 
discussed in Ref. |2(| 

The Curie temperature is estimated within the mean- 
field approximation for a multi-sublattice material by 
solving the system of coupled equations^Sl 




FIG. 1: Schematic representation of the L2i structure 
adapted by the full Heusler alloys. The lattice consists of four 
interpenetrating fee sublattices with the positions (0, 0, 0) and 
(|, i, §) for the Ni and (±, \, \) and (f , f , f ) for the Mn and 
Sn, respectively. 

TABLE I: Lattice parameters, magnetic moments and Curie 
temperatures in NigMnSn at ambient pressure and applied 
pressure of 16 GPa. For comparison, the magnetic moment 
obtained with full potential FLAPW-GGA method is pre- 
sented. 





a=6.022 A[P=0 GPa] 


a=5.821 A[P=16 GPa] 


Ni 


0.21 


0.19 


Mn 


3.73 


3.47 


Sn 


-0.05 


-0.05 


Total 


4.09,4.10 a 


3.81 


T c [calc] 362,373 b 


400 


T c [expt] 360 c , 342 d , 338° 




a Ref(3| 






b Refl28 






c RefH 






d ReflW 






"RefEl 







where (s v ) is the average z component of sj^ and Jq U = 
Sr ^or- Eq.© can be represented in the form of eigen- 
value matrix problem 

(0 - TI)S = (7) 

where M „ = j^—Jq U , I is a unit matrix and S is the 
vector of (s u ). The largest eigenvalue of matrix gives 
the value of Curie temperature^ 

III. RESULTS AND DISCUSSION 

We will subdivide the discussion of the influence of 
the pressure on the electronic properties of Ni2MnSn in 
two parts. First, we present a detailed study of the 
low-pressure region where an experimental information 
is available (We extend this interval up to ~20 GPa). 
In particular we verify the monotonous increase of the 
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Curie temperature with increasing pressure in this region. 
Then we consider a much larger interval of the variation 
of the lattice parameter to study the occurrence of the 
non-monotonous behavior of the Curie temperature. 



A. Low pressure region 

We begin with the discussion of the effect of pressure 
on the electronic structure. In Fig. [2 we compare the 
density of states calculated for the ambient pressure and 
the applied pressure of 16 GPa. As expected, the pressure 
leads to the broadening of the bands that stems from the 
decreasing inter-atomic distances and, therefore, increas- 
ing overlap of the atomic states. One of the consequences 
of the band broadening is the trend to the decrease of the 
magnetic moments. This trend is demonstrated in tabled 
and Fig. |3] In Fig. we present a detailed information 
on the atomic and total magnetic moments for the range 
of pressures up to 20.6 GPa. The dependence of the Mn 
magnetic moment on pressure can be well represented by 
a linear function with a negative slope. The behavior of 
the induced moment of Ni is more peculiar. The depen- 
dence deviates strongly from the straight line and shows 
weak oscillations in the high-pressure part of the curve. 
This behavior reflects the details of the pressure depen- 
dence of the band structure, in particular, the form of the 
DOS in the energy region close to the Fermi level and the 
character of the Mn-Ni hybridization. Since these weak 
oscillations do not play noticeable role in the issues we 
focus on in this paper we do not further investigate their 
origin. The induced moment on Sn has the direction 
opposite to the direction of the Mn moment. Its value 
decreases slowly with increasing pressure. The spin po- 
larization at the Fermi level shows very weak pressure 
dependence. 

Thus the decreasing lattice constant produces a clear 
trend to a monotonous decrease of the atomic magnetic 
moments. For our purpose of the investigation of the 
pressure dependence of the Curie temperature it is impor- 
tant to relate the increasing band width and decreasing 
magnetic moments to the properties of the inter-atomic 
exchange interactions. 

In the spirit of the Heisenberg model of localized mo- 
ments one expects that decreasing atomic moments pro- 
duce the trend to the decrease of the inter-atomic ex- 
change interactions by the factor of /Mq where M p is 
the atomic moment at pressure P and Mo is the moment 
at the ambient pressure. Correspondingly, one expects 
the trend to decreasing Curie temperature resulting from 
decreasing atomic moments. 

An opposite monotonous trend to increasing inter- 
atomic exchange interactions is produced by increasing 
electron hopping and, as a result, more efficient me- 
diation of the exchange interactions between magnetic 
atoms. The competition of two opposite trends opens 
possibility for both increase and decrease of the Curie 
temperature with applied pressure as well as for a non- 
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FIG. 2: Upper panel: Spin projected density of states of 
Ni2MnSn for ambient pressure and applied pressure of 16 
GPa. Lower panel: The same for the case of Mn-Ni atomic 
interchange. 
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FIG. 3: Pressure dependence of magnetic moments in 
Ni 2 MnSn. 



monotonous pressure dependence in a larger pressure in- 
terval. 

In Fig. 0^, we present the calculated inter-atomic ex- 
change parameters of Ni2MnSn for pressures of and 16 
GPa. For comparison, a zero-pressure result of previ- 
ous calculation is also presented. At both pressures the 
patterns of inter-atomic exchange interactions are very 
similar. This similarity involves both Mn-Mn and Mn- 
Ni exchange interactions. The Mn-Mn interactions are 
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long-ranged reaching beyond the 8th nearest neighbor- 
hood distance and have the RKKY-type oscillating char- 
acter. The inter-sublattice Mn-Ni interaction behaves 
very differently. A sizable interaction takes place only 
between nearest neighbors. Note that FigQJi does not 
present all calculated exchange parameters: the exchange 
parameters have been evaluated up to the inter-atomic 
distances of 8.7a that corresponds to about 70 coordi- 
nation spheres. The absolute value of the parameters 
decays quickly with increasing interatomic distance. In 
Fig.|3t>, we demonstrate the convergence of the calculated 
Curie temperature with respect to increasing number of 
the atomic coordination spheres. The main contribution 
to T c comes from the interaction between atoms lying 
closer than 3a. After 5a no sizable contribution is de- 
tected. 

In FigHk, we compare our exchange parameters for 
NiaMnSn at zero pressure with the exchange parame- 
ters calculated recently by Kurtulus et al£* Kurtulus et 
al. used TB-LMTO-ASA method and local spin-density 
approximation (LSD A). The inter-atomic exchange pa- 
rameters were evaluated using the real-space approach 
by Liechtenstein et al~ This approach and the frozen- 
magnon technique employed in the present paper are 
equivalent to each other. In the real-space method by 
Liechtenstein et al. the inter-atomic exchange parame- 
ters are calculated directly whereas in the frozen-magnon 
approach they are obtained by the Fourier transforma- 
tion of the magnon dispersion. 

In Tabled we present the MFA estimation of the Curie 
temperature. It is in good agreement with available ex- 
perimental values overestimating them somewhat. An 
overestimation of the Curie temperature is a usual feature 
of the MFAi 30 i 31 i 32 It arises from the property that the 
MFA expression for T c corresponds to an equal weight- 
ing of the low- and high-energy spin-wave excitations. A 
better weighting of the magnetic excitations is provided 
by the random-phase approximation (RPA). 32 However, 
in the case of the lattices with high atomic coordina- 
tion numbers and in the cases of the magnon disper- 
sion deviating strongly from a simple sinusoidal form the 
both MFA and RPA can give similar values of the Curie 
temperature making the MFA estimation reliable i 3 ^* 3 ^ A 
good agreement of the theoretical and experimental T c 
values shows that the MFA is well applicable for the given 
system. 

The pressure dependence of the interatomic exchange 
parameters is presented in Fig. |SJ The corresponding 
Curie temperature is shown in inset in Fig [HJ The anal- 
ysis shows that the leading contribution into Curie tem- 
perature is given by the Mn-Mn exchange interactions 
within the first three coordinations spheres. The num- 
bers of the atoms in these spheres are 12, 6 and 24, re- 
spectively for the first, second and third spheres. The 
exchange parameters corresponding to the second and 
third coordination spheres increase monotonously with 
increasing pressure determining the increase of the Curie 
temperature (Fig|HJ. Thus, the increase of the exper- 
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FIG. 4: a)Upper panel: Interatomic exchange parameters 
of Ni2MnSn for ambient pressure and applied pressure of 16 
GPa. Lower panel: The same for the case of Mn-Ni atomic in- 
terchange. Zero pressure comparison of both Mn-Mn (+) and 
Mn-Ni (★) exchange parameters with Refl2ij is also shown, 
b) Pressure variation of the Curie temperature with increas- 
ing number of coordination spheres with and without Mn-Ni 
atomic interchange. 



imental Curie temperature with pressure in the corre- 
sponding pressure region^ is well confirmed by the cal- 
culations. In terms of the competition of the two opposite 
trends discussed above this result means a stronger effect 
of the increasing hopping compared with the effect of de- 
creasing atomic moments. Such a behavior is expected 
for large inter-atomic distances. 

Note that in the interval from 9 GPa to 16 GPa we 
obtain a flat feature in the pressure dependence of the 
Curie temperature. This behavior is in a correlation with 
the recent experiment of Kyuji et al~ who obtained the 
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FIG. 5: Pressure dependence of the first four nearest neigh- 
bor Mn-Mn exchange parameters and two Mn-Ni exchange 
parameters. 
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Mn-Mn distance. Upper inset shows the variation of induced 
Ni moment with Mn-Mn distance. Lower inset shows Erot 
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Curie temperature increase from 338K at the ambient 
pressure to 395K at 12 GPa. At ^7 GPa they obtained 
a decrease of the pressure gradient that can be put into 
correspondence to the theoretical flat feature. 

Both measured and calculated Curie temperatures are 
in good agreement with the empirical interaction curve 
for the corresponding region of the Mn-Mn distances 
(Fig. El- 



B. High pressure region 

To verify the non-monotonous pressure dependence 
of T c predicted by the interaction curve we extended 
the calculation to smaller Mn-Mn distances down to 
3.09A The calculated magnetic moments are presented 
in Fig. The Mn moment decreases with the reduction 
of Mn-Mn distance. An interesting feature is obtained 
at dMn-Mn = 3.416A where the value of the magnetic 
moment changes discontinuously. To study the origin 
of the discontinuity we employed the fixed-spin-moment 
methodSLiSiiS that allows the calculation of the total en- 
ergy as a function of the spin moment for a given lattice 
parameter. The corresponding curves are presented in 
the inset in Fig. |S| For large Mn-Mn distance the curve 
has one minimum corresponding to a high-spin state. In 
the region of the discontinuity the curve has two local 
minima revealing the presence of a metastable state. At 



the point of discontinuity the minimum corresponding to 
the low-spin state becomes lower. With further decrease 
of the lattice volume the minimum corresponding to the 
high -spin state disappears. At dMn-Mn — 3.09A the 
magnetic moment vanishes and the ground state of the 
system becomes a Pauli-paramagnet. 

In Fig. [7J the first four nearest neighbor Mn-Mn ex- 
change interactions are presented for the broad interval 
of the Mn-Mn interatomic distance. The pressure region 
discussed in the preceding section corresponds to the last 
three points in the plot. Three of the four leading param- 
eters show nonmonotonous behavior that is reflected in 
the nonmonotonous behavior of the Curie temperature 
(Fig. ISJ). The absolute values of the 2nd, 3rd and 4th 
neighbor Mn-Mn parameters first increase with pressure, 
reach their maxima at the Mn-Mn distances in the region 
from about 4. OA to about 3.6A and decrease strongly 
with further decrease of the Mn-Mn distance. There are 
some weak peculiarities in the behavior of the exchange 
parameters at the region after the discontinues transition 
such as additional local extrema in the 1st, 3th and 4th 
neighbor parameters. They, however, compensate each 
other and the Curie temperature has only one maximum 
at the Mn-Mn distance of about 3.8 A (Fig©. 

The non-monotonous behavior of the exchange param- 
eters and T c can be interpreted as a result of the com- 
petition of two opposite monotonous trends appearing 
with the variation of the Mn-Mn distances discussed in 
the previous section. In the low-pressure region the in- 
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FIG. 7: The first four nearest neighbor Mn-Mn exchange 
interactions in NLMnSn as function of Mn-Mn interatomic 
distance. The ellipses show the region of the metamagnetic 
transition. 

fluence of the increasing hopping prevails while in the 
high-pressure region the influence of decreasing magnetic 
moments becomes more important. 

Qualitatively, the calculated pressure dependence of T c 
in the broad pressure interval is in agreement with the 
Kanomata's empirical interaction curve. Indeed, we ob- 
tained non-monotonous pressure dependence character- 
ized by one maximum separating the regions of positive 
and negative pressure gradients. The low-pressure part 
of the calculated dependence is in reasonable quantita- 
tive agreement with the Kanomata's interaction curve. 
For smaller lattice volumes the calculated Curie temper- 
ature decreases faster than it is prognosticated by the 
interaction curve. The calculations predict the discon- 
tinuity in the pressure dependence of the Curie temper- 
ature of NLMnSn which is absent in the empirical in- 
teraction curve. The extension of the measurements to 
higher pressures is desirable to verify the predictions of 
the calculations. 



C. Atomic inter-sublattice interchange 

The calculations of the Curie temperature of NigMnSn 
discussed in the preceding sections are in good correla- 
tion with measured T c values for the range of pressures 
studied experimentally. A detailed numerical compari- 
son shows, however, that the theoretical pressure deriva- 
tive, dT c /dP, estimated as 3.22 K/GPa is substantially 
smaller than the experimental estimation of 7.44K/GPa 
obtained by Gavriliuk et ali^ To verify the role of the 
atomic interchange between Ni and Mn sublattices ob- 
served experimentally^ we performed calculation for a 



FIG. 8: Schematic representation of Kanomata's empirical 
interaction curve and the calculated Curie temperature as a 
function of the interatomic Mn-Mn distance in NLMnSn. In- 
set shows pressure variation of T c in the low pressure region. 
Small rectangles present the Curie temperatures of the cor- 
responding compounds at ambient pressure. The attached 
arrows show the sign of dT c /dP. (The experimental informa- 
tion is taken from Ref lioh . 

model system where the atoms of the Mn sublattice are 
interchanged with the atoms of one of the Ni sublat- 
tices. Although this model is a strong simplification of 
the experimental situation it allows the investigation of 
the trends resulting from the Mn-Ni interchange. 

With Mn-Ni interchange we obtain a substantial dif- 
ference in the electron structure of the system. The cor- 
responding DOS and magnetic moments for ambient and 
applied pressure of 16 GPa are presented in FigEl an d 
table ITT1 In this case the Mn states hybridize differently 
with the states of two Ni atoms. As a result the magnetic 
moments of the Ni atoms assume different values. The 
total magnetic moment per formula unit decreases from 
3.50/xb at ambient pressure to the 3.09^s at the applied 
pressure of 16 GPa. The decrease of the total magnetic 
moment is mostly the result of the reduction of the Mn 
moment. Note that at the pressure of 16 GPa the rela- 
tive variation of the total magnetic moment is two times 
larger than in Ni2MnSn without Mn-Ni interchange. The 
change in the shape of the 3d peaks and broadening of 
the bands are similar to those for the system without 
swapping. 

At zero pressure the pattern of exchange parameters 
(Fig. and resulting Curie temperature (Table [nj are 
very similar to the case without Mn-Ni interchange. 
However, the situation is different at applied pressure of 
16 GPa. Both Mn-Mn and Mn-Ni 1 nearest-neighbor ex- 
change parameters increase substantially. The remaining 
exchange parameters show small pressure dependencies. 
The interaction of the Mn moment with the moment of 
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TABLE II: Lattice parameters, magnetic moments and Curie 
temperatures at ambient pressure and applied pressure of 16 
GPa for Ni2MnSn with Mn-Ni atomic interchange. 

q=6.022 A[P=0 GPa] a=5.821 A[P=16 GPa] 



Np 021 0.26 

Ni 2 0.08 0.05 

Mn 3.24 2.80 

Sn -0.04 -0.03 

Total 3.50 3.09 

T c 340 562 



the second Ni atom is slightly reduced. The substantial 
increase of the leading exchange parameters with pres- 
sure results in considerable change of the Curie tempera- 
ture from 400 K at ambient pressure to 562 K at 16 GPa. 
Assuming a linear variation of T c with pressure we esti- 
mate the pressure derivative, dT c /dP, as 12.5 K/GPa. 
This value of dT c /dP exceeds strongly the corresponding 
value for the system without Mn-Ni atomic interchange. 
Since the number of the swaped Mn and Ni atoms in 
our model is much larger than in the samples measured 
the calculated dT c /dP cannot be directly compared with 
the experimental pressure derivative. Important, how- 
ever, that the Mn-Ni atomic interchange increases the 
pressure derivative of T c that gives an explanation for 
the low value of the theoretical pressure derivative in the 
case of the system without swapping. A detailed study of 
the influence of the inter-sublattice atomic interchange 
on the electron properties of the Heusler systems is an 
interesting extension of the present study. 



IV. CONCLUSION 



In conclusion, we have systematically studied the pres- 
sure dependence of exchange interactions and Curie 
temperature in full Heusler alloy Ni 2 MnSn within the 
parameter-free density functional theory. We show that 
the character of the pressure dependence of the exchange 
interactions is a consequence of the complex interplay of 
competing trends in the electronic properties of the sys- 
tem. In agreement with experiment, the Curie temper- 
ature increases with increasing pressure in the pressure 
region studied. Extending our theoretical study onto a 
larger pressure interval we obtained non-monotonous T c 
dependence and the presence of the metamagnetic tran- 
sition. The T c behavior in the whole pressure interval 
is in qualitative correlation with the Kanomata's empir- 
ical interaction curve. In the low-pressure region there 
is good quantitative agreement between calculated val- 
ues and the prediction of the empirical rule. The Mn- 
Ni atomic interchange is shown to increase the pressure 
derivative of the Curie temperature that suggests the 
physical mechanism for the improved agreement between 
the experimental and theoretical estimations of this pa- 
rameter. 
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